close all; clear all; clc; 


psi1=0.1417;psi2=1.65;

psi1b=0.1407;psi2b=1.05;

a = load('agrid2.txt');
a = exp(a);

udef = max(0,psi1+psi2*log(a));
udefb= max(0,psi1b+psi2b*log(a));

figure;
plot(a,udef,'r',a,udefb,'b');
